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A  METHOD  OF  DETERMINING  SOME  AERODYNAMIC  COEFFICIENTS 
FROM  SUPERSONIC  FREE  FLIGHT  TESTS  OF  A  ROLLING  MISSILE 


ABSTRACT 

Based  on  the  assumption  that  the  aerodynamic  forces  acting  at 
any  point  on  a  lifting  surface  are  linearly  dependent  upon  the  local  angle 
of  attack  at  that  point,  the  theory  of  pure  rolling  motion  is  applied  to  the 
experimentally  observed  motion  of  a  finned  missile  in  order  to 

1.  determine  how  well  the  actijal  motion  is  represented  by  the 
theory, 

2.  determine  the  suitability  of  the  Aerodynamic  Range  to  the  pro¬ 
posed  free -flight  roll  technique,  and 

3.  determine  certain  aerodynamic  coefficients  associated  with  the 
motion. 

Values  for  the  aerodynamic  coefficients  are  also  derived  using  the 
linearized  supersonic  theory  and  are  compared  with  those  obtained  from 
experiment. 

The  results  indicate  that  the  actual  motion  is  well  represented  by 
the  theory  to  within  the  small  experimental  errors  and  that  excellent  re¬ 
producibility  of  the  aerod3mamic  coefficients  in  roll  is  obtained.  Further  - 
more  the  results  show  a  fairly  good  correlation  with  the  linearized  theory 
considering  the  degree  of  approximation  associated  with  this  theory  when 
applied  to  aerodjmamic  surfaces  of  16%thickness  as  employed  on  the  mis¬ 
siles  used  in  the  reported  tests. 


4 


BALLISTIC  RESEARCH  LABORATORIES 


SYMBOLS  AND  COEFFICIENTS 


a 

A,B 

b 

B 


c 


radius  of  body  of  missile 
area  of  fin  or  wing  panel  (see  Figure  1) 
arbitrary  constants  depending  on  boundary  conditions 
wing  span 

|/m^  -  1 

wing  chord 
damping  constant 


I 


drag  coefficient 

roll  moment  derivative  due  to  canted  surface 


6  in  rad. 


roll  moment  derivative  due  to  rolling  velocity, 
roll  moment  coefficient,  • 

qApb 

diameter  of  body 


ec„ 


a  (2^) 


slope  of  the  lift  curve  at  any  point  x,y  on  fin 


K. 


K. 


R 


I 

M 


m 


LjV 

Lp/pV 


rolling  moment  due  to  angle  of  cant 

rolling  moment  due  to  rolling  velocity,  damping  moment 

2 

axial  moment  of  inertia  (slugs  -  ft  ) 

Mach  number 
mass  (slugs) 
number  of  canted  fins 
total  number  of  fins 
rolling  velocity  (^,  rad/sec 
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P.E 

Eb 

2V 

2V 

5 
s 
V 
x,y 

z 

a 

6 
/S 

.r 

0 

S 


steady  state  rolling  velocity 
probable  error 

fin  tip  helix  angle 

steady-state  fin  tip  helix  angle 


steady  state  rolling  velocity  (°/ft) 

linear  velocity  of  missile  along  trajectory,  ft/sec 

Cartesian  coordinates,  y  along  fin  span,  x  along  fin  chord 

linear  distance  along  the  free-flight  range 

angle  of  attack 

angle  of  cant  on  one  fin 

tip  mach  cone  angle 

polar  coordinates  shown  in  figure  3 

free  stream  air  density,  — 

f^ 

angle  of  roll 
s  -  0' 


Superscripts 

"dot"  over  symbol  is  time  derivative  of  variable 

."prime"  mark  above  sjrmbol  is  the  derivative  with  respect 
to  Z,  distance  along  range. 
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INTRODUCTION 

A  continually  guided  supersonic  missile  is  possible  only  after  roll  control  has  been  obtained,  and  the 

difficulties  in  achieving  this  are  well  known.  The  design  of  a  reliable  roll  stabilization  system  requires  an 

accurate  knowledge  of  the  dynamic  as  well  as  the  static  aerodynamic  coefficients  associated  with  the  rolling 

motion.  Since  the  Aerodynamic  Range  appeared  to  offer  promise  as  a  new  tool  for  the  determination  of  both 

static  and  dynamic  coefficients  associated  with  aerodynamic  configurations,  an  attempt  to  determine  these 

critical  rolling  coefficients  seemed  to  be  a  valuable  initial  undertaking.  Some  of  the  ballistic  contributions 

to  aerodynamics,  as  determined  in  the  Aerodynamic  Range  are  presented  by  Charters^,  and  a  detailed  de- 

2 

scription  of  the  Range  is  given  by  Charters  and  Thomas  . 

In  this  preliminary  report  the  theory  of  pure  rolling  motion  is  applied  to  the  experimentally  observed 
motion  in  order  to  obtain  the  aerodynamic  coefficients  associated  with  the  motion.  The  fit  of  the  equation  of 
rolling  motion  to  the  range  data  is  obtained,  and  a  comparison  is  made  of  the  experimentally  obtained  aero¬ 
dynamic  coefficients  to  those  predicted  by  a  linearized  theory. 

It  is  not  intended  that  this  paper  present  the  results  of  a  finished  research  program,  but  rather  to 
discuss  the  basic  theory  of  motion,  the  general  test  procedure,  the  method  of  data  reduction,  some  specific 
results  from  the  experimental  technique,  and  the  experimental  accuracy  that  is  attained. 


GENERAL  EXPERIMENTAL  TECHNIQUE 

A  suitable  aerodynamic  test  vehicle  is  gun-laimched  down  the  Aerodynamic  Free -Flight  Range.  The 
special  firing  technique  employed,  results  in  very  small  pitching  and  yawing  displacements  of  the  test  vehicle 
during  the  subsequent  flight.  This  vehicle  may  be  a  model  of  any  aerodynamically  stable  missile  configura¬ 
tion  having  one  or  more  pairs  of  fins  set  at  a  differential  angle  of  incidence  (cant)  to  induce  a  rolling  velocity. 
During  flight,  as  it  passes  each  of  twenty -five  stations  accurately  located  to  0.001  feet  along  a  218 -foot  dis¬ 
tance,  the  missile  is  photographed  (shadowgraphs)  simultaneously  in  the  vertical  and  horizontal  planes  by 
means  of  an  extremely  short -duration  spark  discharge.  The  time  of  occurrance  of  the  discharge  for  9  of  the 
stations  is  measured  to  within  10  ^  seconds.  From  the  photographs,  the  spatial  orientation  of  the  missile  is 
obtained,  from  which  the  linear  displacements  (to  0.001  ft.  accuracy)  and  angular  displacements  in  roll  (to 
+.7  degree  accuracy)  are  determined  at  the  25  discrete  distances  along  the  trajectory,  as  well  as  at  9  discrete 
values  of  time. 

With  these  data  obtained  from  the  Aerod3mamic  Range  and  the  analyses  in  the  following  sections,  the 
aerodynamic  coefficients  associated  with  the  rolling  motion  may  be  determined  and  compared  with  the  theo¬ 
retical  values  calculated  on  the  basis  of  linearized  supersonic  theory. 


1  -  "Some  Ballistic  Contributions  to  Aerodynamics",  A.  C.  Charters,  Journal  of  the  Aeronautical  Sciences, 

Vol.  14,  No.  3,  p.  155,  March  1947. 

2  -  'The  Aerodynamic  Performance  of  Small  Spheres  from  Subsonic  to  High  Supersonic  Velocities",  A.  C. 

Charters  and  R.  N.  Thomas,  Journal  of  the  Aeronautical  Sciences,  Vol.  12,  No.  4,  p.  468,  October  1945, 
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THEORY 

A.  Equation  of  Motion  of  a  Missile  in  Pure  Roll. 

If  a  S3nninetrical  winged  and/or  finned  missile  with  canted  surfaces  is  in  free  flight  and  has  an  angu¬ 
lar  motion  of  pure  roll  (  yawing  and  pitching  displacements  negligible),  the  missile  may  be  considered  as 
being  acted  upon  by  two  roll  moments.  One  moment,  tending  to  increase  the  roll  velocity,  results  from  the 
lift -force  couple  produced  by  the  canted  surfaces  in  linear  forward  motion  and  the  second  moment,  tending  to 
decrease  the  roll  velocity,  is  the  damping  moment  produced  by  the  lift -force  couples  induced  on  all  the  aero¬ 
dynamic  surfaces  by  the  roll  velocity  itself. 

1.  Roll  Moment  of  a  Missile  Due  to  Canted  Fins. 

The  roll  moment  induced  by  n ^  differentially  offset  (canted)  wing  or  tail  panels,  may  be  expressed  as  the 
integrated  product  of  the  aerod3mamic  lift  force  and  the  moment  arm  of  this  force  about  the  missile  axis.  A 
double  integration  is  convenient  because  a  finite  wing  moving  through  air  at  supersonic  velocities  experiences 
a  variable  pressure  over  the  surface  within  the  region  bounded  by  the  wing-tip  Mach  cone.  Consequently, 
from  Figure  1  the  roll  moment  may  be  written  as: 

^6 

where  /  \ 

n  =  number  of  identical,  canted  wing 

or  tail  panels  on  the  missile 

j  =  angle  of  cant 

Then 

a)  L 

6 

b)  L.  =  K  (2) 

®  *8 

2.  Roll  Moment  Due  to  Rolling  Velocity  -  Damping  Moment.  When  a  missile  rolls  about  its  own  axis, 
the  lifting  surfaces  are  subjected  to  a  normal  velocity  vector  induced  by  the  rolling  velocity.  The  induced 
velocity  results  in  an  induced  angle  of  attack  on  the  surface  which  varies  over  the  surface  as  a  linear  function 
of  the  radial  distance  from  the  rotational  axis.  If  the  roll  damping  moment,  neglecting  the  viscous  damping 
on  the  fuselage,  ^  is  defined  as  the  roll  moment  induced  by  and  opposing  or  damping  the  rolling  velocity,  then 


The  effect  of  viscous  damping  of  a  missile  fuselage  may  be  estimated  from  the  performance  of  spinning 
bodies  of  revolution  as  determined  in  the  Aerodynamic  Range.  The  numerical  value  of  the  viscous  damp¬ 
ing  coefficient  of  our  missile  is  approximately  0.16%  of  the  aerodynamic  damping  coefficient  obtained  due 
to  the  fins.  Consequently  the  effect  of  body  damping  may  then  be  neglected  for  the  purposes  of  this  report. 


nj  PV 


b/ 


ydxdy 


c  b/ 


f  f  PV 

“  I  1  (  - ^  — TT- 

x=0  y=a 


ydxdy 


(1) 


=  slope  of  the  lift  curve  at  any  point  x,y 
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FIGURE  2;  SCHEMATIC  SHOWING  ANGLE  OF  ATTACK  OF 
ANY  ELEMENT  OF  A  ROLLING  SURFACE  FOR 
PURE  ROLLING  MOTION. 
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this  aerodynamic  damping  moment  (with  reference  to  Figures  1  and  2)  may  be  expressed  by  a  double  integra¬ 
tion  as  follows: 


where  p^  is  the  induced  angle  of  attack  on  the  lifting  surfaces,  and  n  is  the  number  of  identical  wing  or 


a) 


the  missile  fuselage. 

c 

2 

=  pV  Pn  1 

f  {«L) 

2 

x=0 

J  '  da  ' 

y=a 

y  d  X  d  y 


(4) 


x,y 


b)  L  =  PV 

^  TV 


(4) 


3.  Dynamical  Equation  of  Roll.  Applied  to  a  symmetrical  missile  in  pure  rolling  flight,  subject  to 
small  pitching  and  yawing  displacements,  the  dynamical  equation  in  roll  becomes  simply; 


10 


L 

P 


(5) 


Here  the  rolling  moments  exerted  by  the  wing  and  tail  surfaces  arising  from  the  pitching  and  yawing  angular 
displacements  and  angular  velocities  have  been  neglected  because  (1)  for  small  pitching  and  yawing  displace¬ 
ments  (mean  square  yaw  less  than  1.5“)  these  moments  are  extremely  small  and  of  second  order  compared 
with  Lj  and  (see  Appendix  A)  and  (2)  they  are  oscillatory  with,  respect  to  their  direction  of  action  and 
as  such  any  small  effect  that  they  may  otherwise  exert  on  the  resulting  rolling  motion  tends  to  cancel  out, 

B,  Solution  of  Dynamical  Equation  of  Roll. 

For  the  missiles  traveling  the  length  of  the  range  between  stations  (218  feet),  the  decrease  in  velo¬ 
city  is  less  than  7  percent  and  for  values  of  flight  Mach  number  in  excess  of  1.6  the  variation  of  dC  /da 

L  ^ 

the  slope  of  the  lift  curve,  resulting  from  this  small  deceleration  is  less  than  about  3.5  percent  from  the 
average  value.  Therefore  little  error  is  introduced  if  the  values  of  and  in  equations  (2)  and  (4) 

are  considered  as  constant  for  any  given  flight.  ® 

Insertion  of  equations  (2)  and  (4)  into  (5)  gives: 


10 


«y 


d  X  d  y 


(6) 


where 


Pn 

~2 


(V) 


10 
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c  b/2 


C  =  ™  j  j  /d.dy 

'  d  « 'x,y 


and  0  =  ^  =  P  rolling  velocity. 
^  dt 


The  approximate  linear  equation  of  motion  of  the  missile  along  the  trajectory  may  be  written: 


where 


mV  =-K_V 

irv 


is  considered  constant  for  the  small  variation  in  velocity  V. 

Equations  (6)  and  (9)  are  now  transformed  to  the  independent  variable,  z,  the  distance  along  the 
trajectory,  ^  by: 

^  =  V  (10) 

dt  dt  dz  dz 

_v2  A  .vi!-  -3- 

2  ,2  ''  dz  dz 

dt  dz 

If  a  prime  is  used  to  denote  a  derivative  with  respect  to  z,  then  equations  (6)  and  (9)  become: 


v* 

V  m 


Substitution  of  (12)  into  (11)  gives: 


"  */ 

0  +0(  A  - 


m  /  I 


a  linear  differential  equation  with  constant  coefficients. 


In  the  Range,  the  distance  along  the  trajectory  may  be  assumed  equal  to  the  distance  along  the  range 
with  negligible  error,  and  the  vector  velocity 

V  =  ^ 

^  dt  dt  dt 


where 


becomes 


dx  dz 

dT“‘^Tt«^ 


V  =  negligible  error, 

dt 
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Defining 


Equation  (12)  becomes: 


Ci=  _!p  .  ^ 


s  = 


0  +C^0  -C^=0 


The  general  solution  of  Equation  (14)  maybe  written  down  immediately  as: 


where 


-C  z 

0=B  +  sz  +  Ae  1 


s  =  ^2  = 

C, 


steady -state  rolling  velocity. 


A  and  B  are  arbitrary  constants  which  depend  upon  the  boundary  conditions  written  as  follows: 
For  z  =  0,  0  =  0  and  0'  =  0' 

z  =  oo  ,  0'  =  0' 


Then  equation  (15)  becomes; 

(0*00  -00^ 

0  =  00  +  0*00  ^  - 


e^l^ 


Equation  (15)  then  should  represent  the  form  of  the  curve  connecting  the  set  of  experimental  data  points  of 
roll  angle  versus  distance  along  the  trajectory  obtained  from  the  firings  in  the  range.  A  curve  of  the  form 
given  in  equation  (15)  then  is  fitted  to  the  data  points  by  successive  applications  of  the  method  of  Differential 
Corrections  and  the  constants  C.  and  C„  as  well  as  K  and  K.  are  determined  experimentally.  Thus 

the  damping  and  rolling  moments  can  be  experimentally  determined  as  a  function  of  missile  design  and  Mach 
number. 

B.  Theoretical  Determination  of  the  Roll-Moment  Coefficients. 

The  values  for  K  and  K  defined  in  equations  (7)  and  (8)  are  determined  experimentally  in  the 

Aerod3mamic  Range  by  the  technique  discussed  in  this  paper.  However,  approximate  values  for  these  two 
coefficients  may  also  be  determined  theoretically  from  the  application  of  the  existing  linearized  supersonic 
theory. 


The  theory  as  applied  here  assumes  the  two-dimensional  valve  of 


in  the  areas  of  the  fins  out¬ 


side  the  tip  region  for  the  calculation  of  K  and  K  .  This  assumption,  of  course,  cannot  give  the  cor- 
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rect  lift  distribution  on  the  wing  within  the  Mach  cone  originating  from  the  intersection  of  the  wing  leading 
edge  with  the  body.  Since  this  region  in  which  the  lift  distribution  is  incorrectly  calculated  is  close  to  the 
body,  the  resultant  error  in  the  rolling  moments  should  be  small. 

The  effects  of  body -wing  interference  when  the  body  is  at  an  angle  of  yaMr  have  also  been  neglected. 
However,  in  the  roll  tests  in  the  Aerodynamic  Range,  it  should  be  noted  that  the  body  yaw  during  flight  down 
the  range  is  always  intentionally  restricted  by  the  launching  technique  to  a  very  small  angle  (mean  square  ydiw 
less  than  1.5')  in  order  for  equation  (6)  to  be  very  nearly  exact.  Any  upwash  on  the  wing  due  to  body -wing 
interaction  resulting  from  the  yawed  body  would  result  in  (1)  small  symmetrical  forces  (lift  forces  in  the 
same  direction  on  both  halves  of  the  fin)  when  the  fins  are  oriented  at  0  or  45  degrees  to  the  plane  of  yaw, 
and  (2)  possibly  very  small  unsymmetrical  forces  due  to  flow  dissimilarities  when  the  missile  is  oriented  in 
roll  to  any  angle  other  than  those  stated  above  referred  again  to  the  plane  of  yaw.  Since  a  roll  moment  can 
obviously  arise  only  from  the  unsymmetrical  forces,  it  may  be  concluded  from  the  qualitative  observations 
above  that  any  roll  moments  resulting  from  body  upwash  phenomena  would  not  only  be  very  small  moments 
but  would  be  rapidly  oscillating  both  in  sign  as  well  as  in  magnitude.  Such  moments  would  contribute  negli¬ 
gible  effect  to  the  rolling  motion  of  the  missile  resulting  from  the  relatively  large  and  continuously  directed 

moments,  L  and  L ,  . 

P  6 

1.  Determination  of  Kj  from  the  Linearized  Supersonic  Theory. 

P  1 

From  the  linearized  two-dimensional  theory  of  Ackeret  ,  the  value  of  the  lift  slope  at  all  points 
on  the  surface  of  an  infinite -aspect-ratio  lifting  surface  is  given  by  the  expression; 


which  is  a  function  of  Mach  number  only.  However,  for  a  finite  wing  or  tail  surface,  the  pressure  distri¬ 


bution  over  the  lifting  surface  included  within  the  tip  Mach  cone  is  influenced  by  the  wing-tip  boundary 

conditions  as  well  as  by  the  variable  angle  of  attack  ( helix  angle)  induced  over  the  surface  by  the 

rolling  velocity.  This  pressure  distribution  is  evaluated  on  the  basis  of  the  three-dimensional  linearized 
2  3  fdC 

supersonic  theory  ,  using  Eward’s  method.  The  value  then  fori _ L  1  for  any  incremental  area 

1.  da  yx,  y 

on  a  rectangular  wing  within  the  Mach  cone  boundary  is. 


a  function  of  position  (x,  y)  on  the  wing  as  well  as  of  the  flight  Mach  number. 

^Liepmann,  H.  W.,  and  Puckett,  A.  E.,  Aerodynamics  of  a  Compressible  Fluid,  Chapter  9,  John  Wiley  and 
Sons,  Incl.,  1947. 


9 

Harmon,  S.  M.  "Stability  Derivatives  of  Thin  Rectangular  Wings  at  Supersonic  Speeds",  NACA  TN  No.  1706, 
November,  1948. 


3 

Eward,  J.  C.,  'Distribution  of  Wave  Drag  and  Lift  in  the  Vicinity  of  Wing  Tips  at  Supersonic  Speeds", 
NACA  TN  No.  1382,  1947. 
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From  Figure  3, 


For  Area  I 


=  (^Jpjarea  I  +(^*p  )  H  jarea  III 


From  equations  (8)  and  (19)  using  dS  in  place  of  dxdy  for  the  differential  area 


I  s 
where  from  Figure  3, 

dS  =  rdiUdr=  x  dx  d  fi 

cos  n  cos  fl 

Inserting  (22)  into  (21)  gives; 


y^dS 


tan 

r 

/K  ]  - 

4Pn  J 

i  v)' 

o 

II 

a. 

ffl 

b 

I 

From  Figure  3, 

and 

tan 

II 

4  Pn  f 

TTB  J 

Ml) 


/ _  (.3in-Vf(!^ - 


y  =  ^  -  X  tan  iu 


lU 


(20) 


(21) 


(22) 


c 

ff^(^  ~  ^  sin'^KB  tan  jU  -  -  x  tan  -  tan^  dfi  d: 

ai  cos^  f/  J 


FIGURE  3-  SKETCH  INDICATING  METHOD  OF 

INTEGRATION  FOR  ROLL  MOMENTS. 
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2.  Determination  of  K  from  Linearized  Supersonic  Theory. 

'a 

The  expression  for  K  is  determined  by  an  integration  of  Equation  (7).  Here  again  Equation  (18) 

% 

gives  the  value  of  the  lift  slope  on  the  wing  or  tail  surface  at  all  points  outside  of  the  area  bounded  by  the 
Mach  cone.  The  lift  slope  for  the  area  within  the  Mach  cone  is  given  by: 
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derived  by  Snow  and  Bonney.  This  value  differs  from  Equation  (19)  because  the  lift  forces  that  give  rise 
to  the  roll  moment  are  a  result  of  a  uniform  angle  of  attack  at  all  points  on  the  surface --namely  the 
angle  of  cant. 

Referring  again  to  Figure  3, 


K. 


area  I  area  II  '''  ^  area  III 

n  6  O' 


(28) 


For  Area  I: 


From  Equations  (7),  (22)  and  (27); 


0 


Pn,  6 
6 


-1  (-) 

c  tan  ^bV 

-  /  / 


jU=0 


/B(|-y) 


] 


J2L 


2 

cos  fi 


d  jU  d  X 


and 


d 


"“a*  (c\ 

b2  I  4 


4.8  B 


(29) 


For  Area  II: 


From  Equations  (18)  and  (7) 


(\)ll 

Q 


Pn,  fi 
0 


B(  I  -  y)  I 

f*  ^  /*  o 


f  f 


x=0 


b  c 

y=2-B 


3  y  dx  dy 


and 


(K  >1,  =  2P.j«( 


For  Area  III: 


2,  3 

c  b  c 


4B  3B 


From  Equations  (18)  and  (7) 


Pn  6 
0 


J 


b  c 
2  ■  B 


3  cy  dy 


(30) 


and 


Pn  fi 
6 

B 


( 


b  c 


2.  3 

c  b  +  c 


2 

a  c 


B  B 


) 


(31) 


^Snow,  R.  M.,  and  Bonney,  E.  A.,  "Aerod5mamic  Characteristics  of  Wings  at  Supersonic  Speeds",  Bumblebee 
Repcfrt  No.  55,  The  Johns  Hopkins  University  Applied  Physics  Laboratory,  March  1947. 
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And  adding  Equations  (29),  (30),  and  (31) 


(K,  ) 

6 


Pn.  6 
0 

B 


(t 


2  c^b 

a  c  - 


8B 


2k 

c  b 
4B  J 


(32) 


3.  Determination  of  the  Conventional  Aerodynamic  Coefficients  in  Roll. 

It  is  desirable  to  express  any  experimentally  or  theoretically  determined  coefficients  in 
terms  of  the  conventionally  defined  aerodynamic  coefficients  of  roll,  namely: 


ec 


s  5 


roll  moment  derivative  due  to 
canted  surfaces 


,  roll  moment  derivative  due  to 
rolling  velocity 


By  definition: 


L.  =  Co 


q  n  8  b  =  Kp  V 
P  0  * , 


and  L  =  C 


«p(  2v) 


qAbn=Ko  Vp 


Where 

Then 


area  of  tail  or  wing  panel  (defined  as  half  of  one  exposed 
tail  or  wing  surface,  see  Fig.  1) 


2K„ 


p  b  A  n  6 
P  a 


(33) 


4K. 


pb  A  n 
P 


These  coefficients  may  also  be  expressed  in  terms  of  the  experimental  constants  and  Cg  through 


use  of  equations  (13a),  (13b),  and  (33)  as, 

s.  = 


2  I  '^2 


P  b  A  n  « 

pa 


(34) 


41 


pb  A  n 
P 


(-1  *  >) 


Comparisons  of  the  experimental  and  theoretical  values  of  Cj  and  Cj  offer  a  means  of  eval- 

p  a 

uating  the  accuracy  of  the  linearized  supersonic  theory  as  applied  here  when  used  to  calculate  these  quantities. 
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It  will  be  noted  from  equation  (32)  that  accurate  measurements  of  the  actual  angle  of  incidence  of  each  tail 

and  wing  panel  is  required  for  the  theoretical  determination  of  Cj;  for  any  missile  tested. 

6 

METHOD  OF  DATA  REDUCTION 


The  experimental  values  of  the  aerod3mamic  coefficients  in  roll  for  the  missiles  tested  are  determined 
by  the  following  general  procedure  for  reduction  of  the  data  obtained  from  the  Aerodynamic  Range. 

A)  Accurate  measurements  of  the  angular  orientation  in  roll  (0)  and  the  linear  position  (z)  of  the 
missile  in  space  are  obtained  from  each  of  the  photographic  plates  taken  during  the  flight  of  the  missile  down 
the  range.  Very  accurate  measurements  of  the  time  at  which  9  of  the  photographs  were  taken  is  determined 
from  the  chronograph  record  and  megacycle  counters. 

B)  The  form  of  the  solution  0  =  0(z)  equation  (15)  of  the  theoretical  roll  equation  is 'fitted"  to  the 
measured  values  of  0  and  z,  and  the  optimum  values  of  the  constants  and  Cg  in  equation  (15)  are 
determined  by  the  method  of  Differential  Corrections. 

Kr 

C)  The  value  of  the  retardation  constant  — fi—  for  the  missile  is  determined  by  existing  techniques 
using  the  measured  values  of  z  and  the  time  t. 

D)  Values  of  the  aerodynamic  coefficients  in  roll  are  determined  and  their  statistical  accuracy 
(Probable  Error)  is  computed. 

The  details  of  these  four  general  steps  will  now  be  discussed. 

A.  Measurements  from.the  Photographic  Plates. 

1.  Angular  Orientation  in  Roll.  In  order  to  determine  the  roll  angle  from  the  shadowgraphs  of  the 
missile  (see  Figure  4)  it  is  first  necessary  to  identify  the  quadrant  in  which  a  particular  reference  fin  or 
wing  occurs.  Therefore  one  particular  fin  must  be  recognizable  in  both  the  vertical  and  horizontal  photo¬ 
graphs  taken  at  each  station.  Identification  of  the  fins  was  obtained,  after  considerable  experimention,  by 
locating  a  small  pin  in  the  trailing  edge  of  each  fin  near  the  tip.  The  shape  of  the  four  pins  differed  as 
seen  in  Figure  4,  and  it  was  always  possible  to  distinguish  at  least  one  of  the  shapes  in  both  photographs 
and  thereby  definitely  fix  the  quadrant  of  roll  orientation.  The  roll  angle  of  the  missile  at  any  station  may 
be  calculated  from  the  coordinates  of  each  of  the  identification  pins  (see  Eq.  35)  as  obtained  from  the 
horizontal  and  vertical  spark  photographic  plates  using  the  usual  methods  of  reduction  of  range  data  for 
obtaining  the  spatial  position  of  a  point  from  its  images  projected  on  the  horizontal  and  vertical  plates. 

(See  Fig.  5). 


<p  =  tan 


AC  -  CD 
M) 


<p  =  tan 


BD  -  AB 
AD 


cj>  = 


tan 


EG  -  GH 
EH 


(35) 


(j)  =  tan 


HF  -  EF 
EH 


•^\  t} 


f /*  ic-^ 


^'jK‘ V\  ^  ”*' 
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FIGURE  5 


20 


BALLISTIC  RESEARCH  LABORATORIES 


The  angular  accuracy  of  the  roll  angle  at  a  station  is  estimated  to  be  ^  .7°  from  a  comparison  of  the 
four  individual  determinations  made  according  to  Eq.  35. 

2.  Linear  Position  of  Missile.  The  location  of  a  reference  point  on  each  station  along  the  range  is 
known  to  .001  feet.  This  reference  point  appears  on  the  photographs  and  the  location  of  the  center  of  gravity 
of  the  missile  in  space  at  each  station  is  accurately  and  simply  determined  by  the  technique  presented  by 
Charters  and  Thomas  \ 

B.  Determination  of  the  Constants  in  the  Roll  Equation. 

The  constants  and  Cg  and  at  the  same  time  the  constants  A,  B,  and  s  are  determined  by  'fitting" 
the  theoretical  form  of  the  roll  equation  (15),  0  =  0(z),  to  the  0,  z  data  obtained  from  a  particular  firing. 
The  process  of  "fitting"  consists  of  first  finding  the  initial  or  approximate  values  of  the  constants  in  equation 
(15)  (designated  as  A^,  B^,  S^,  and  )  and  second  applying  the  method  of  Differential  Corrections  to  deter 
mine  the  optimum  values  of  the  constants. 

1.  Method  of  Differential  Corrections. 

Since  the  roll  equation  contains  an  exponential  term,  no  method  of  least  squares  has  been  developed 
to  handle  this  type  of  equation  directly  (at  least  not  to  the  knowledge  of  the  authors).  However,  the  roll 
Equation  (15)  may  be  expanded  by  a  Taylor  Series  about  the  initial  values  of  the  constants, 


-C,  z  "*^1  ^  ‘^1  ^ 

0  =  B^  +  s^z  +  A^e  o  +AB  +  zAs+e  o  AA  -  zA^e  o  " 

AC^AA  '^1^  AaAA  2,.  ■‘^1  ^  .AC.,^ 

zeo  1  -zeo  1  +zAeo(l)  +... 

2  2  °  2 

in  which  all  terms  in  AC^  and  A  A  of  order  2  or  greater  may  be  neglected  if  AC^  and  A  A  are  small. 

-C  z  -C,  z 

0  -  0  =  A  B  +  z-A  s  +  e  o  AA-Aze  oAc, 

m  0  o  1 


where 


0  =(B  +  sz+Ae  o);  initial  value  of  roll  angle 

0  0  0  0  ’ 


0  =  measured  value  of  roll  angle 

m  ^ 

is  a  linear  relation  and  the  differential  corrections  (  AB,  As,  A  A,  and  AC^)  to  the  initial  value  of  the  con¬ 
stants  may  be  obtained  by  the  usual  least  squares  process.  The  corrected  constants 


*1  ■  = 

Bj  =  B  .  iB  Cj  -  Cj  .  iCj 

o 

may  in  turn  be  used  as  the  initial  values  and  the  process  repeated  until  the  sum  of  the  squares  of  the 


^Charters,  A.  C.  and  Thomas,  R.  N.,  'The  Aerodynamic  Performance  of  Small  Spheres  from  Subsonic 
to  High  Supersonic  Velocities",  JOURNAL  OF  THE  AERONAUTICAL  SCIENCES,  Vol.  12,  No.  4,  P.  468, 
October  1945. 
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1  2 

residuals'^  approaches  a  minimum.  If  the  process  is  not  convergent,  then  either  the  roll  equation  does 

not  represent  the  motion  with  sufficient  exactness,  or  the  initial  values  for  the  constants  are  not  well  deter¬ 
mined. 

Therefore,  the  successful  application  of  the  method  of  Differential  Corrections  depends  upon  a  fairly 
accurate  determination  of  the  initial  values  of  the  constants.  The  techniques  employed  to  obtain  them  are 
discussed  in  the  following  section. 

2.  Determination  of  the  Initial  Values  for  the  Roll  Constants. 

(a)  Steady  State  Rollii^  Velocity,  (s^).  Differentiation  of  the  roll  equation  (15), 

I  “C,  z 

0  =  s  -  Cj  A  e  (36) 

shows  that  as  z  increases,  0',  the  roll  velocity,  approached  a  constant  value,  s.  The  approximate  value  of 
0'oo  (designated  s^)  may  be  determined  from  the  range  data  in  the  following  manner. 

O 

The  average  rolling  velocity"^  0  *  between  successive,  closely  spaced  stations  is  plotted  for  the  value 

of  z  midway  between  the  corresponding  stations.  The  values  of,  0'^,  0^, - 0^  taken  from  the  resulting 

plot  at  equal  but  arbitrary  intervals  of  distance  (Az)  along  the  range  are  then  substituted  into  eqmtion  (36) 

-Cl  A 
1  Az 

0'  =  s  -  Cj  A  e 

-C^2Az 

0|  =  s  -  A  e 


-C 


^n-l=  "  -  A  e 


0^  =  s  -  A  e 


1  (n-1)  A  z 


-C 


1  n  A  z 


where  the  ratio  of  the  differences  between  successive  values  of  0'  becomes; 

40' 

n 


A  0' 


n-1 


0'  -0'  , 

n  n-1 

'  -  0' 

n-1  ^  n-2 


-C^  Az 

=  e  =  constant 


^Residual  here  is  defined  as  the  difference  between  the  roU  angle  of  the  missile  as  measured  and  the 
angle  calculated  from  the  roU  equation  at  any  distance,  z,  along  the  range  where  the  missile  was  photo¬ 
graphed. 

o 

^If  the  sum  of  the  squares  of  the  residuals  decreases  with  successive  applications  of  the  method  of  Dif¬ 
ferential  Corrections,  the  process  is  convergent. 

3 

This  is  determined  simply  by  dividing  the  0  between  adjacent  stations  by  the  corresponding  distance  between 
the  stations. 
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or 


A0' 


n-1 


-C^Az 

e 


(37) 


From  Equation  (37)  a  table  of  A  0'  for  any  interval  may  be  calculated  since  any  difference  ^  0'j^  can  be  found 
from  the  previous  known  difference  ^  i‘  Then  the  roll  rate  is  given  by  ^ 

0\-0\  +  4^2  +4^'+.  ...  40^ 


and  the  approximate  steady  state  rolling  velocity  by: 

oo 


s 

o 


^0 

n 

n  =  2 


t 


In  practice,  the  interval  of  distance,  A  z,  used  in  reading  from  the  curve  of  0  versus  Z  was 

t 

50  feet,  and  that  used  in  predicting  4  0  at  distances  greater  than  the  length  of  the  range  was  100  feet. 
Values  of  4  0  were  predicted  until  they  became  of  negligible  magnitude. 


(b)  Damping  Constant,  (C^  ).  An  initial  value  for  the  damping  constant  of  the  roll  equation  and  a 

0 

check  on  the  value  of  s^  as  previously  determined,  maybe  obtained  frgm  Equation  (36)  rewritten  in  terms 
of  the  approximate  constants  as; 


Log  S  =  log  (  C,  A^j  -  z  =  Const  -  z  (38) 

where 

S2=s^-0 

A  plot  of  log  S  against  Z,  using  the  values  of  s^  and  0  previously  calculated  from  the  range  data, 
should  yield  a  straight  line  with  slope,  ~  the  damping  constant.  If  the  resultant  plot  is 

not  quite  a  straight  line,  the  value  of  s^  may  be  slightly  dhanged  until  a  constant  slope  for  the  plot  is  ob¬ 
tained. 

(c)  Boundary  Constants,  and  Although  initial  values  for  A  and  B  maybe  obtained  directly 
from  equations  (15)  and  (38)  for  z  =  0,  it  was  found  to  be  abetter  procedure  to  solve  for  them  from  equation 
(15)  alone  by  a  least  squares  process  using  the  range  data  and  the  approximated  values  for  s^  and  . 

3.  "Fit"  of  the  Roll  Equation.  The  final  equation  representing  the  rolling  motion  of  the  missile 
is  equation  (15)  written  with  the  optimum  values  of  the  constants  derived  from  the  data  by  successive 
applications  of  the  method  of  Differential  Corrections  until  the  sum  of  the  squares  of  the  residuals  ap¬ 
proaches  a  minimum.  The  fit  of  the  final  roll  equation  to  the  range  data  may  then  be  evaluated  by  a  com¬ 
parison  of  the  statistical  magnitude  of  the  final  residuals  with  the  magnitude  of  the  estimated  error  in  the 


1 


This  method  was  suggested  by  Dr.  E.  J.  McShane,  Professor  of  Mathematics,  University  of  Virginia. 
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measurement  of  the  roll  angle  from  the  photographic  plates.  The  statistical  magnitude  of  the  final  resi¬ 
duals  is  given  by  the  usual  statistical  quantity,  Probable  Error 


where 


P.E. 


(39) 


A0  -  residual 
n  =  number  of  observations 

If  the  probable  Error  of  the  residuals  is  of  the  same  order  of  magnitude  as  the  estimated  error 
in  measurement  of  the  roll  angles  from  the  shadowgraphs  (+  .7“),  then  the  final  roll  equation 
represents  the  best  possible  equation  that  can  be  written  to  fit  the  data  and  as  such  represents  the  true 
rolling  motion  of  the  missile  within  experimental  accuracy. 


C. 


Determination  of  the  Value  of  the  Retardation  Constant 


The  retardation  constant  defined  by  equation  (9b)  is  proportional  to  the  coefficient  of  drag  of  the  mis- 

2 

sile.  The  constant  then  is  determined  in  the  conventional  manner  described  in  where  it  is  assumed  that  the 
time  is  given  by  a  power  series  in  the  distance  measured  from  the  center  of  the  range  (z^)  as: 

2 

t  =  a  +  a,  (z  -  z  )  +  ao  (z  -  z  )  + . 

0  1  ^  o'  A'  o' 

The  coefficients  a  ,  a  ,  a_,  etc.,  are  determined  from  the  experimental  measurements  by  the  method  of 

0  J.  ^ 

Least  Squares  and  the  drag  coefficient  is  in  turn  determined  from  these  coefficients. 

D.  Aerodynamic  Coefficients  and  the  Statistical  Precision 

The  values  of  the  constants  C.  and  of  the  roll  equation,  as  determined  by  successive  applica¬ 
tions  of  the  method  of  Differential  Corrections,  may  be  considered  as  the  best  values  that  the  experimental 

data  can  yield.  The  roll  coefficients  Cj  and  C„  are  then  determined  simply  by  equation  (34). 

P  ^6 

The  effectiveness  parameter,  p^b  ,  the  steady-state  tip  Helix  angle,  maybe  obtained  from  the 

W 


linearized  theory  by 


sb 


w 

I 

\.  p  J 


(radians) 


2V  2  (57.3) 

and  from  the  experimental  data  by  substitution  of  equations  (13a)  and  (13b)  into  equation  (40),  whence; 


(40) 


2V 


r 


2(57.3) 


mj 


^Scarborough,  J.  B.,  "Numerical  Mathematical  Analysis",  The  Johns  Hopkins  Press,  1930. 

2 

Charters,  A.  C.  and  Thomas,  R.  N.,  'The  Aerodynamic  Performance  of  Small  Spheres  from  Subsonic  to 
High  Supersonic  Velocities",  JOURNAL  OF  THE  AERONAUTICAL  SCIENCES,  Vol.  12,  No.  4,  P.  468 
October  1945.  ’ 
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and  using  equation  (16)  gives: 


P  b 
s 

2V 


sb 


r 


2  (57.3) 


1 


1 

1  +  — 
CiHly 


(41) 


The  accuracy  of  the  final  aerodynamic  coefficients  is  again  expressed  by  the  statistical  quantity, 
Probable  Error,  (P.E.),  dependent  upon  the  P.E.  of  the  coefficients  in  the  roll  equation.  The  general  equations 
for  the  probable  error  of  any  function  of  independent  quantities  whose  P.E.’s  are  known  is  given  by:^  , 


R 


where 


(42) 


and  R,  r^  r^. 


r^  are  the  P.E.’s  of  Q,  q^,  qg,  —  q^  respectively,  determined  by  Equation  (39). 


EXPERIMENTAL  TESTS 

The  purpose  of  the  experimental  program  is  not  to  present  a  finished  research  program  but  rather 
to  investigate  the  suitability  of  the  Aerodynamic  Range  to  the  free  -flight  roll  technique  presented  in  this 
paper.  Specifically,  the  purpose  may  be  restated  in  three  parts: 

1)  To  determine  the  accuracy  of  agreement  of  the  theoretical  roll  equation  with  the  experimental 
data  measured  in  the  Aerodynamic  Range  as  well  as  to  determine  the  accuracy  of  the  reduction  methods 
used  to  obtain  the  roll  constants  (aerodjmamic  coefficients)  from  the  data. 

2)  To  determine  the  ability  of  the  combined  Range  and  technique  to  accurately  reproduce  aerodynamic 
coefficients. 

3)  To  compare  the  experimental  and  theoretical  values  of  the  aerodynamic  coefficients. 

A.  Experimental  Apparatus  and  Procedure 

1.  Aerodynamic  Range 

The  Aerodynamic  Range  and  the  test  procedure  are  described  briefly  in  the  introduction  to  this 

2 

paper.  Further  details  are  presented  by  Charters  and  Thomas  and  Figure  6.  A  unique  modification  to 
the  equipment  is  the  special  "railed"  gun  (shown  in  Figure  7)  used  for  launching  winged  and/or  finned  mis¬ 
siles.  In  the  operation  of  the  gun,  the  projectile  is  guided  by  its  body  in  the  central  cylindrical  portion  of 
the  barrel,  the  fins  or  wings  riding  freely  in  the  slots. 

2.  Design  of  Missiles  for  Firing  Program 

To  fulfill  the  purposes  of  the  program,  a  single  missile  design  was  chosen  consisting  of  the 
simplest  configuration  from  a  construction  standpoint  as  well  as  from  a  standpoint  of  amenability  to 

^Scarborough,  J.  B.,  'Numerical  Mathematical  Analysis",  The  Johns  Hopkins  Press,  1930 

^Charters,  A.  C.  and  Thomas,  R.  N.,  "The  Aerodynamic  Performance  of  Small  Spheres  from  Subsonic  to 
High  Supersonic  Velocities",  JOURNAL  OF  THE  AERONAUTICAL  SCIENCES,  Vol.  12,  No.  4,  P.  468, 

October  1945. 
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Figure  6.  Aerodynamic  Range. 
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Figure  7.  Special  Railed  Gun. 
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theoretical  supersonic  aerod3mamic  calculations.  The  missile  design  is  shown  in  Figure  8,  and  consists  of 
a  cone -cylinder  fuselage  with  four  rectangular  fins  in  cruciform  arrangement.  The  fin  surfaces  are  wedge 
shaped  with  16%  maximum  thickness  and  were  constructed  as  shown  in  the  fin  detail  of  Figure  8,  for  reasons 
of  simplicity.  On  each  model,  two  of  the  tail  fins  were  canted. 

For  two  of  the  models,  the  fuselage  was  of  aluminum  with  wings  of  Tobin  bronze.  Bronze  wings 
were  used  to  obtain  (1)  a  higher  moment  of  inertia  in  roll  in  order  to  decrease  the  rate  of  acceleration  in 
roll  and  correspondingly  lengthen  the  curved  portion  of  the  measured  roll  displacement  -  distance  curve, 
and  (2)  a  higher  mass  in  order  to  decrease  the  linear  deceleration  of  the  missile  and  result  in  more  nearly 
constant -velocity  flight.  The  remaining  two  models  were  machined  entirely  from  Tobin  bronze  to  further 
increase  the  moment  of  inertia  in  roll.  A  summary  of  the  missile  design  is  given  below. 


Missile 

No. 

Average 

Angle  of 

Cant  on 
two  Tail  fins, 

Material 

Moment  of 
Inertia  in 
roll  -  I 
(slugs  -ft  j 

Mass 

of 

Missile  - 

m 

(slugs) 

1 

2.029° 

Aluminum 
fuselage 
+  Bronze 
fins 

1.240  X  lO'^ 

1.224  X  10“' 

2 

4.288° 

Aluminum 
fuselage 
+  Bronze 
fins 

1.235 

1.221 

3 

2.108° 

All  bronze 

2.186 

3.068 

4 

4.454° 

All  bronze 

2.210 

3.067 

B.  Results  and  Discussion 

The  four  missile  models  (Table  I)  were  gun  launched  at  an  intended  Mach  number  of  about  1.7. 
The  actual  Mach  number  that  resulted  in  the  fo;^  flights  both  at  the  beginning  and  end  of  the  trajectory  as 
well  as  the  air  density  and  velocity  of  somd  in  the  range  are  recorded  in  Table  II. 

TABLE  II 


Missile 

No. 

Mach  No. 
at  begin¬ 
ning  of 
traj  ectory 

Mach  No. 
at  end  of 
trajectory 

Mach  No. 
at  center 
of  timing 
data 

Vel.  of 
Sound 
in  Range 
(Ft/Sec) 

Air  Density 
in  Range  „ 
(Slugs/Ft  ) 

1 

1.790 

1.675 

1.746 

1135 

2.297  X  10"° 

2 

1.752 

1.635 

1.694 

1131 

2.296 

3 

1.621 

1.575 

1.592 

1121 

2.357 

4 

1.665 

1.616 

1.637 

1120 

2.357 
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1,  Rolling  Motion 

The  angular  orientation  in  roll  of  the  models  as  measured  from  the  shadowgraphs  taken  at  each 
spark  station  along  the  trajectory  is  shown  by  the  points  in  Figure  9.  In  this  same  figure  are  plotted  the 
equations  of  the  curves  which  best  fit  the  data  as  obtained  by  the  methods  discussed  in  the  previous  section. 

The  difference  between  the  curves  for  missiles  number  2  and  1  as  well  as  number  4  and  3  are  a 
result  of  the  differences  in  the  angles  of  cant  employed.  The  differences  between  the  curves  for  missiles 
number  2  and  4  as  well  as  number  1  and  3  are  due  primarily  to  the  differences  in  the  roll  moment  of  in¬ 
ertia  (see  Table  I). 

The  use  of  a  greater  moment  of  inertia,  it  was  felt,  would  result  in  a  better  determination  of  the 
roll  coefficients  because  of  the  slower  roll  acceleration  rate  and  linear  deceleration  rate.  However^  no 
such  variation  in  experimental  accuracy  with  moment  of  inertia  can  be  noted  in  the  results  and  the  ac¬ 
curacy  appears  dependent  only  upon  the  exactness  of  the  measurements  obtained  from  the  spark  photo¬ 
graphs. 

The  excellent  "fit"  of  the  roll  equation  to  the  experimental  points  which  cannot  be  observed  in 
Figure  9  because  of  size  limitations  is  given  in  Table  HI  as  the  statistical  quantity.  Probable  Error: 

TABLE  in 

Model  12  3  4 

P.  E.  .26°  .73°  .65°  .57° 

Since  the  Probable  Error  is  small  and  of  the  same  order  of  magnitude  as  the  estimated  error  in  roU  angle 
measuLrement,  it  may  be  concluded  that  the  roll  equation  represents  the  actual  motion  of  the  models  tested 
to  within  the  experimental  error.  It  therefore  appears  that  the  rolling  motion  is  unaffected  by  the  pitching 
motion  for  the  mean  square  yaws  encountered  in  the  tests  (see  Table  IV). 

The  optimum  values  of  the  constants  in  the  roll  equation  were  obtained  by  the  methods  outlined 
in  the  Data  Reduction  Section  and  are  given  with  their  Probable  Error  in  Table  IV  as  determined  from  the 
reduction. 

TABLE  IV 


Model 

1 

2 

3 

4 

A 

660° 

1386 

882 

2285 

B 

-621° 

-1298 

-838 

-2079 

^1 

.0166 

.0169 

.0096 

.0098 

P.E. 

.24% 

.34 

2.34 

.53 

s 

14.90 

32.02 

15.57 

32.64 

P.E. 

.04% 

.05 

.57 

.15 

Cl  s 

.2463 

.5418 

.1502 

.3211 

P.E. 

.24% 

.34 

2.41 

.55 

(^) 

.  2.7  x  lO"^ 

2.8 

1.2 

1.2 

P.E. 

.7% 

.1 

1.1 

.4 

Yaw 

1—^ 

CO 

o 

.34 

.35 

1.15 

No.  of 
obser. 

14 

21 

16 

18 

30 
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2.  Aerodynamic  Coefficients 

The  experimental  values  of  aerodynamic  coefficients  for  the  models  tested  may  be  obtained  from 
the  data  in  Table  IV  and  Equation  34.  The  theoretical  values  of  the  coefficients  for  the  models  are  obtained 
from  Equations  26,  32  and  33.  The  values  for  the  coefficients  are  given  in  Table  V  and  their  very  good 
quality  is  again  indicated  by  the  statistical  factor,  Probable  Error,  the  significance  of  which  is  discussed 
in  the  Meiliod  of  Reduction  Section. 

The  experimental  values  of  Cj  and  (given  in  Table  V)  and  the  theoretical  curves  are 

P  ^ 


plotted  respectively  in  Figures  10  and  11  against  the  average  flight  Mach  number.  (See  Table  II).  It  is 
seen  from  the  figures  that  the  values  for  model  3  do  not  fall  in  line  with  those  for  the  other  models.  It  may 
be  seen  from  Table  V  that  the  P.E.’s  for  model  3  are  more  than  four  times  as  large  as  the  largest  P.E.  for 
the  other  three  models.  This  is  a  positive  indication  that  the  reduction  of  model  3  is  weak,  in  comparison. 

It  is  further  noted  in  Figure  9,  that  model  3  is  lacking  stations  at  the  beginning  of  the  range  where  data  are 
required  for  a  good  determination  of  C^.  The  dashed  line  in  the  figures  represents  the  theoretical  slope 
at  an  average  Mach  number  drawn  thru  the  weighted  average  of  the  coefficients  at  this  average  Mach  number 
(See  "Internal  Consistency  and  Reproducibility  of  the  Data,  and  Linearity"  Sect.  4  following). 

The  weighted  average  experimental  value  of  Cj  at  the  reduced  Mach  number  as  indicated  by  the 

P 


dashed  curve  in  Figure  10  is  21.0%  greater  than  the  value  calculated  from  the  linearized  theory  as  applied 
in  this  paper  (no  body -wing  interference)  and  the  corresponding  value  of  Cp  in  Figure  10  is  only  7.6% 


greater  than  the  theoretical.  (See  Table  VI) 

Concerning  this  deviation  of  the  experimental  values  from  the  theoretical  it  will  be  remembered  that 
the  airfoil  section  used  was  of  triangular  shape  (single  wedge),  with  about  16% maximum  thickness  at  the 
trailing  edge.  A  comparison  of  the  theoretical  lift  of  an  infinite -aspect -ratio  wing  of  such  cross  section  ob¬ 
tained  by  the  linearized  theory  with  that  obtained  from  the  exact  theory  using  the  oblique  shock  equations 
results  in  approximately  a  15%  lower  value  for  the  linearized  lift  at  a  Mach  number  in  the  neighborhood  of 


1.7. 

In  view  of  the  above  results  one  is  tempted  to  speculate.  The  difference  between  the  experimental  and 

the  theoretical  values  for  Cj  might  immediately  be  attributed  to  the  fact  that  the  fins  were  of  such  large  thick 

P 

ness  that  the  linearized  theory  simply  failed  to  predict  the  correct  magnitude  of  pressures  acting  upon  them. 

However,  from  Figure  10,  the  fact  that  the  experimental  values  of  Cj  deviate  only  7.6%  from  the  calcula- 

6 

ted  values  (instead  of  about  21.0%  as  was  the  case  for  C„  )  must  obviously  modify  this  first  observation  in 

P  /dC  V 

the  following  manner.  Let  us  assinne  that  at  any  point  on  the  fin  the  actual  values  of  the  lift  slope  ^  j 

for  the  two  roll  moments  are  equal  except  for  a  small  deviation  in  the  tip  region  (see  Figure  12b)  as  is  in¬ 
dicated  by  the  theory  when  the  rectangular  fins  are  considered  attached  to  an  infinite  wall,  and  let  us  con- 

sider  the  theoretical  distributions  shown  by  the  solid  lines  in  Fig.  12b.  Then,  if  the  above  as- 

sumption  is  reasonably  correct,  since  depends  upon  the  value  of  y  in  the  integration  (see  Figure  1) 

P 

and  Cj  upon  y,  the  respective  deviations  of  the  experimental  results  of  Cj  and  from  the  calculated 

5  P  5 
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values  indicate  that,  in  addition  to  the  magnitude,  the  actual  distribution  of  pressure  on  the  fin  surface 
differs  appreciably  from  the  assumed  theoretical  distribution;  and  it  may  be  inferred  that  this  distribution 
of  pressure  must  be  somewhat  distorted  from  the  calculated  pattern  in  a  manner  qualitatively  shown  in 
Figure  12a  with  greater  pressures  acting  over  the  area  in  the  region  near  the  wing  tips.  Such  a  pressure 
distribution  would  result  in  a  somewhat  lopsided  lift-slope  distribution  over  the  fin  as  indicated  in  Figure 
12b  and  may  be  a  combined  result  of  the  large  fin  thickness  employed,  (particularly  critical  at  the  fin  tip) 

and  the  wing-body  interference. 

3.  Effectiveness  Parameter 

The  experimental  values  for  the  effectiveness  parameter  (obtained  from  equation  (40),  Table  IV, 
and  the  physical  measurements)  and  the  theoretical  values  for  this  parameter  (obtained  from  the  equations 
(26),  (32),  and  (41))  are  given  in  Table  V  together  with  their  Probable  Error  as  obtained  from  the  reduction. 
Figure  13  gives  the  theoretical  cimve  for  the  effectiveness  parameter  per  6  vs  Mach  number  together 
with  the  experimental  points.  It  is  noted  that  the  experimental  points  are  approximately  11.3%  lower 

than  the  theory.  See  Table  VI.  This  difference  is  due  to  the  ratio  of  Cjj  to  Cj  (see  Eq.  41);  these 

5  P 

coefficients  were  discussed  in  paragraph  2.  (Aerodynamic  Coefficients). 


TABLE  V 


Model 

1 

2 

3 

4 

Mach  No. 

1.746 

1.693 

1.592 

1.637 

(th.) 

P 

.449 

.460 

.483 

.472 

(ex.) 

P 

.548 

.559 

.546 

.564 

P.  E.* 

.24% 

.33 

2.32 

.52 

(th.) 

6 

.716 

.740 

.794 

.769 

(ex.) 

8 

.778 

.806 

.784 

.802 

P.  E.  * 

.24% 

.34 

2.41 

.55 

^s^  (th.) 

2V 

.0283 

.0602 

.0303 

.0633 

Ps^  (ex.) 

2V 

.0247 

.0531 

.0260 

.0546 

P.  E.  * 

.04% 

.05 

.59 

.15 

♦The  low  values  for  the  Probable  Errors  indicate 
that  the  coefficients  are  rather  well  determined  from 
the  reductions. 


Internal  consistency  or  Reproducibility  of  the  Data,and  Linearity 

The  internal  consistancy  or  reproducibility  of  the  data  was  estimated  by  reducing  the  experimental 
p  b  / 

C.  ,  Cj  ,  and  — to  a  common  Mach  number  of  M  =  1.665  by  the  theoretical  value 

p  ^8 


4. 

values  of 
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(i2b) 


FIG.  12  -  DISTRIBUTION  ON  THE  FIN  PANEL  OF  (a)  THE 
COEFFICIENTS  OF  PRESSURE  DUE  TO  CANT 
AND  ROLLING  VELOCITY  AND  (b)  OF  THE 
LIFT  SLOPE. 
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of  the  slopes  at  M  =  1.665.  These  reduced  values  are  given  in  Table  VI. 
The  weighted  average  was  obtained  from  ^ 


and  the  Probable  Error  of  the  weighted  average  from 


Weighted  Average 
1 


0.6745 


where 


W 


n  =4 


V  =  Diff.  of  coeff.  at  same  Mach  No.  from  average. 


These  values  of  Probable  Error  are  a  measure  of  the  internal  consistancy  or  reproducibility  of  the 

data  and  their  small  values  indicate  that  the  technique  is  quite  promising. 

The  theoretical  values  of  the  coefficients  at  the  common  Mach  number  are  also  given  in  Table  VI 

together  with  their  deviation  from  the  experimental  weighted  average. 

In  plotting  C,  vs  M  it  was  assumed  that  Cn  was  linear  with  6  ;  this  assumption  may  now  be 
6 

investigated  by  plotting  Cjj  vs  6  where 


It  is  seen  in  Figure  14  that  the  assumption  of  linearity  is  consistant  with  the  experimental  results. 

TABLE  VI 


Model 

Mach  No 

Ps^  / 

P 

5 

2V 

1 

1.665 

.566 

.820 

.01241 

2 

1.665 

.565 

.821 

.01247 

3 

1.665 

.530 

.746 

.01215 

4 

1.665 

.557 

.787 

.01218 

Weighted  Average 

.564 

.816 

.01242 

P.E.  of  Wt.  Av. 

.27% 

.55% 

.26% 

Theory 

.466 

.758 

.014 

Deviation  of  Weighted  Average 

from  Theory 

21.0% 

7.6% 

-11.3% 

In  Figure  15  the  effectiveness  parameter  is  plotted  vs  6  and  it  is  seen  to  be  linear,  in  good  agreement 
with  theory. 

^Scarborough,  J.  B.,  "Numerical  Mathematical  Analysis",  The  Johns  Hopkins  Press,  1930 
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CONCLUSIONS 

1.  The  tests  establish  the  validity  of  the  roll  eqxiation  for  accurately  representing  the  true  motion 
of  the  missile  as  well  as  the  accuracy  of  the  reduction  methods  employed  since  the  probable  error  of  the 
residuals  of  the  "fitted"  equation  is  a  maximum  of  .73°  and  is  of  the  same  magnitude  as  the  estimated  error 
in  the  roll  angle  measurements. 

2.  The  tests  show  that  the  roUing  moment  coefficient  due  to  cant  and  the  effectiveness  parameter 
depend  linearly  on  the  angle  of  cant  for  angles  of  cant  less  than  4.5°. 

3.  The  tests  show  that  the  experimental  values  of  the  coefficients,  and  ,  as  obtained  from  a 

P  d 

16%  thick  wing  are  respectively  21.0%  and  7.6%  larger  than  the  values  calculated  from  the  three-dimen- 

Pgt)/ 

sional  linearized  theory,  and  that  the  experimental  value  of  6  is  11.3%  less  than  theory. 

4.  The  tests  establish  the  excellent  reproducibility  of  the  aerodynamic  coefficients  as  shown  by  the 

weighted  Probable  Errors  of  .27%  for  ,  .55%  for  Cj  ,  and  .26%  for  Yyj  ^  when  the  values  of  the 

P  5 

coefficients  are  reduced  to  a  common  Mach  number. 
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APPENDIX  A 

A.  Rolling  Moments  Due  to  Angle  of  Attack  a  and  Angle  of  Yaw  ^ 

These  two  moments  are  identical  in  nature  because  of  the  symmetry  of  the  missile  surfaces. 
The  magnitude  of  such  a  moment  is  given  by; 


S  ,  b 
w,  t 


where  S  is  the  area  of  the  wing,  tail,  or  both.  Then 
w,t 


^  (■ 


\  I  PV^ 

■j  V  2 


where  the  expression  for 


has  been  obtained  from  Reference  4  and  p  now  is  replaced  by  the 


ratio  v/V,  where  v  is  the  sideslip  velocity.  In  eqmtion  (I),  "A"  is  defined  as  the  aspect  ratio  of  the 
surface  considered.  Now  when  a  ,  the  angle  of  attack,  is  replaced  by  its  equivalent  w/V,  where  w  '  is  the 
component  of  missile  velocity  in  the  vertical  plane;  equation  I  becomes: 


r 

2 

r.  2  1 

1 

/  1  -  _ 

3  +  \ 

_ 

'  AB  3A^B^  ' 

PS  ,b 
w,t 


(vw)  (II) 


But  (vw)  is  the  product  of  two  very  small  disturbance  velocities  compared  with  V  and  as  such  the  moment 
is  of  second  order  and  negligible  for  a  rolling  missile. 

B.  Rolling  Moments  Due  to  Pitching  and  Yawing  Angular  Velocities 

Inspection  immediately  shows  that  these  coefficients  also  depend  upon  the  product  of  two  very 
small  quantities  (qv)  and  (rw)  respectively.  Therefore  such  rolling  moments  are  negligible  and  need  no 
further  consideration. 
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APPENDIX  B 

REDUCTION  OF  ROLLING  MOTION  OF  MISSILE  4 

1.  DATA 

The  shadowgraphic  plates  obtained  from  the  operative  stations  in  this  firing  were  placed  on  an 
illuminated  grid  consisting  of  20  divisions  per  inch  where  measurements  were  made  of  (1)  the  locations  of 
the  identification  pins  (See  Figure  5,  a,  b,  c,  d, . etc.)  and  (2)  the  linear  position  of  the  missile  with  re¬ 

spect  to  the  datum  point  in  the  photograph.  From  the  measurements  of  the  pins,  the  coordinates  of  the  pins 
may  be  determined  from  the  geometry  at  the  station  and  thus  the  roll  orientation  may  be  calculated  from 
Equation  35.  These  data  of  roll  angle  and  linear  position  are  given  in  Table  VII  and  are  plotted  in  Figure  9. 

TABLE  VTI 


Sta. 

Z  (ft) 

0° 

40° 

Z  (ft) 

0 

Z  ft. 

o/ft. 

5 

0 

205.76 

107.41 

9.6238 

11.16 

4.8119 

6 

9.6238 

318.17 

146.99 

11.0271 

13.33 

15.1374 

7 

20.6509 

460.16 

167.04 

11.0347 

15.14 

26.1682 

8 

31.6856 

627.20 

169.34 

9.9845 

16.96 

36.6778 

9 

41.6701 

796.54 

144.75 

7.8427 

18.46 

45.5914 

10 

49,5128 

941.29 

157.40 

8.1573 

19.30 

53.5914 

11 

57.6701 

1098.69 

329.47 

15.8524 

20.78 

65.5963 

13 

73.5225 

1428.16 

271.14 

12.0822 

22.44 

79.5636 

14 

85.6047 

1699.30 

281.53 

11.9903 

23.48 

91.5999 

15 

97.5950 

1980.83 

295.97 

12.0661 

24.53 

103.6280 

16 

109.6611 

2276.80 

125.34 

5.0372 

24.88 

112.1797 

17 

114.6983 

2402.14 

127.84 

4.9560 

25.79 

117.1763 

18 

119.6543 

2529.98 

128.02 

4.8814 

26.23 

122,0905 

19 

124.5357 

2658.00 

130.25 

5.0552 

25.76 

127.0633 

20 

129.5909 

2788.25 

985.49 

36.0178 

27.36 

147,5998 

21 

165.6087 

3773.74 

687.80 

23.9589 

28.71 

177.5882 

24 

189.5676 

4461.54 

240.10 

8.1727 

29.38 

193.6540 

25 

197.7403 

4701.64 

2.  Initial  values  for  the  constants. 

From  the  differences  in  the  roll  angles  and  the  linear  distances  between  successive  stations  ap¬ 
proximate  values  of  0'  (deg/ft)  at  a  point  midway  between  the  station  may  be  obtained.  These  values  of 
0'  and  Z  are  also  given  in  Table  VII  and  are  plotted  in  Figure  16. 
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From  the  data  given  in  Table  VII,  an  initial  value  of  s  may  be  obtained  by  the  method  indicated  in 
the  section  on  "Method  of  Data  Reduction"  2,  (a)  as  illustrated  in  Table  VIIL 

TABLE  VIII 


z 

0' 

A  0' 

0 

10.00 

10.00 

14.40 

100 

24.40 

5.20 

200 

29.60 

1.859 

300 

.671 

400 

.242 

500 

.087 

600 

.031 

700 

.011 

800 

.004 

900 

2^  =  32.505 


Ratio 


2.769 


s 

0 


An  initial  value  of  C  may  be  obtained  from  the  slope  of  the  logarithmic  plot  of  (s^  -  0 ' )  vs  z 
(see  section  on  "Method  of  Data  Reduction"  2,  (b),  Table  IX  and  Figure  17. 


TABLE  IX 


and  from  Figure  17, 


Z 

0' 

%  - 

0 

10.00 

22.505 

50 

18.60 

13.905 

100 

24.40 

8.105 

150 

27.60 

4.905 

200 

29.60 

2.905 

4  (In  s  -  0') 

A  Z 

In  22.505  -  In  2.905 

,11373  -  1.06643 

200 

200 

.010236 
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For  Z 


0  from  equation  (15),  an  initial  value  of  A  may  be  obtained. 


(s.  -  0' 


A  = 


22.505 

.010236 


2198.6° 
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When  Z  = 
B 

0 


0,  0 

=  ^0 


is  equal  to  B+A  from  which  an  initial  value  for  B  is  obtained. 


A 


o 


=  205.76  -  2198.61 
=  -  1992.84 

These  initial  values  for  the  constants  are  then  substituted  into  equation  (15)  and  values  for  0  are 
determined,  which  when  compared  with  the  measured  values  yield  the  residuals  at  all  stations.  (The  5th 
range  station  at  Z  =  20  ft.  was  taken  as  the  origin). 

By  the  Process  of  Differential  Corrections,  (see  section  "Method  of  Data  Reduction"  B,  1.)  new  values 
for  the  constants  may  be  obtained. 

TABLE  X 


Sta. 

AB 

A  s 

A  A 

-c  z 

-c.  z 

Z  « 

e  ^ 

-zAe 

d0 

xlO 

xlO"^ 

5 

1 

0 

1 

0 

.0100 

6 

1 

.096238 

.90618 

-1.917385 

.8520 

7 

1 

.208509 

.80947 

-3.675258 

2.0337 

8 

1 

.316856 

.72301 

-5.036797 

.4826 

9 

1 

.416701 

.65276 

-5.980344 

-  .2713 

10 

1 

.495128 

.60242 

-6.557904 

.2298 

11 

1 

.576700 

.55415 

-7.026280 

-1.3931 

13 

1 

.735225 

.47115 

-7.616011 

-4.7240 

14 

1 

.858047 

.41634 

-7.835991 

-5.8101 

15 

1 

.975950 

.36825 

-7.901664 

-8.2936 

16 

1 

1.096611 

.32547 

-7.847147 

-10.4756 

17 

1 

1.146983 

.30911 

-7.795038 

-12.9006 

18 

1 

1.196543 

.29382 

-7.729616 

-12.5386 

19 

1 

1.245357 

.27950 

-7.652862 

-11.7044 

20 

1 

1.295909 

.26541 

-7.562058 

-14.7953 

21 

1 

1.656087 

.18357 

-6.683948 

-20.1296 

24 

1 

1.895676 

.14365 

-5.987121 

-23.3452 

25 

1 

1.977403 

.13212 

-5.743968 

-23.5488 

18 

16.185923 

8.436380 

-110.549392 

-146.322100 

20.556000 

5.028149 

-112.251725 

-215.882537 

5.134326 

-44.089211 

-34.171729 

763.887375 

1036.121366 

.899218 

6.001327 

-2.557996 

-12.843722 

-84.307071 

.468688 

1.180298 

7.723962 

34.407683 

-6.141633 

84.933581 

137.484728 

-  .426238 

.089981 

2.249480 

-1.527194 

-2.140147 

57.446128 

-42.964797 

24.999500 

1.210253 

-4.785710 

A  B  -81.934704 
A  s  12.390916  x  10 
A  A  81.883242 
A  C^-  3.954305  x  10 
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TABLE  XI 


Sta. 

AB 

As 

AA 

ACi 

z 

-c  z 
e 

-CiZ 

-zAe 

d0 

-2 

10 

10 

5 

1 

0 

1 

0 

.0400 

6 

1 

.096238 

.90964 

-1.996384 

-  .4999 

7 

1 

.206509 

.81610 

-3.826605 

.0039 

8 

1 

.316856 

.73211 

-5.290130 

-1.4690 

9 

1 

.416701 

.66360 

-6.306075 

-1.6768 

10 

1 

.495128 

.61431 

-6.936386 

.4210 

11 

1 

.576700 

.56691 

-7.455766 

-1.0870 

13 

1 

.735225 

.48504 

-8.132535 

-2.1645 

14 

1 

.856047 

.43066 

-8.407373 

-1.2416 

15 

1 

.975950 

.38275 

-8.518207 

-1.6392 

16 

1 

1.096611 

.33988 

-8.499753 

-1.6550 

17 

1 

1.146983 

.32344 

-8.460166 

-3.1825 

18 

1 

1.196543 

.30804 

-8.405501 

-1.9323 

19 

1 

1.245357 

.29359 

-8.338028 

.2344 

20 

1 

1.295909 

.27935 

-8.255652 

-2.4564 

21 

1 

1.656087 

.19598 

-7.401556 

-2.0667 

24 

1 

1.895676 . 

.15482 

-6.692978 

-2.1567 

25 

1 

1.977403 

.14285 

-6.441746 

-1.4262 

18 

16.185923 

8.639050 

-119.364841 

-25.265300 

20.556000 

5.235000 

-122.357940 

-27.884490 

5.297827 

-  48.685516 

-  9.571296 

892.690541 

190.743340 

.899218 

6.001327 

-2.533389 

-15.022926 

-5.165477 

.479947 

1.151541 

8.603281 

2.554709 

-6.631380 

101.136922 

23.199535 

-.422138 

.082101 

2.261533 

.374165 

-2.503267 

63.530527 

10.268967 

27.545742 

1.234922 

-.037686 

-5.403122 

A  s  1.341582  X  lO""^ 

AA  5.397988 
AC^  .030517  x  10'^ 

The  new  values  now  become 
B  +  A  B  =  -2074° 

s  +  A  s  =  32.63 

A  +  A  A  =  2280° 

+  A  =  .009841 

The  above  process  is  repeated  until  the  corrections  become  negligibly  small,  and  the  sum  of  the  residuals 
squared  reaches  a  minimum.  (See  Tables  XI  and  XII.  ) 

In  the  second  application  of  the  Differential  Corrections,  the  corrections  in  and  s  were  .suf¬ 
ficiently  small  to  make  a  third  application  unnecessary. 
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The  final  values  are: 

B  = 

2079° 

s  = 

32.64  (°/ft.) 

A  = 

2285° 

o 

II 

.009838 

11 

u 

.321112 

TABLE  Xn 

Sta 

Residuals 

Residuals 

Residuals 

from  Original  from  First 

from  Second 

Constants 

Correction 

Correction 

5 

.0100 

.0400 

-  .0800 

6 

.8520 

-  .4999 

-  .2066 

7 

2.0337 

.0039 

.5492 

8 

.4826 

-1.4690 

-  .7373 

9 

-  .2713 

-1.6768 

-  1.5083 

10 

.2298 

-  .4210 

.5966 

11 

-  1.3931 

-1.0870 

.0042 

13 

-  4.7240 

-2.1645 

-  .8776 

14 

-  5.8101 

-1.2416 

.1561 

15 

-  8.2936 

-1.6392 

-  .1587 

16 

-10,4756 

-1.6550 

-  .1152 

17 

-12.9006 

-3.1825 

-  1.6246 

18 

-12.5386 

-1.9323 

-  .3598 

19 

-11.7044 

-  .2344 

1.3491 

20 

-14.7953 

-2.4564 

-  .8418 

21 

-20.1296 

-2.0667 

-  .4749 

24 

-23.3452 

-2.1567 

-  .5983 

25  2 

-23.5488 

-1.4262 

.0730 

X  = 

=  2426.04 

49.1402 

10.1249 

The  Probable  Error  is  given  by  equation  (39) 

P.E.,^  = 

o 

.6745 

-  =  .5736° 

The  Probable  Errors  of  and  s  can  be  obtained  from  the  equations: 


P.E.„  = 

^1 

P  E 

0 

^c 

P.E. 

P.E. 

s 

0 

W 

s 

Where  W  and  W 
c,  s 


can  be  computed  from  the  last  Differential  Correction  Least  Squares. 


*^1 

=  .00052  = 

.53% 

P.E. 

s 

=  .0501  = 

.15% 

^2  = 

C^  s 

P.E.p 

^2 

=  ^  (%P.E, 

■cf  * 
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APPENDIX  C 

VARIATION  OF  THE  AERODYNAMIC  DERIVATIVES  WITH  MACH  NUMBER 

The  rolling  motion  of  a  rectangular  fin  stabilized  missile  at  a  Mach  number  of  1.7  and  the  associated 

1 

aerodynamic  coefficients  namely  Cj^  ,  Cj^  and  -^j  6  are  given  in  a  previous  paper.  Recently  an  initial 
study  of  the  variation  of  these  coefficients  with  Mach  number  has  been  undertaken  for  the  same  aerodynamic 
configuration  (See  Tables  I  and  E)  with  the  firing  of  three  roimds  near  M  =  2. 5.  Unfortunately  only  one  of 
the  three  roimds  has  excellent  statistical  quality;  however  the  results  from  the  range  firings  and  subsequent 
reduction  for  all  three  models  are  given  herein. 

As  indicated  in  a  previous  paper  ^  the  equation  expressing  the  rolling  motion  is 

-C.Z 

0  =  B  +  SZ  +  Ae  (1) 


where 

0  =  roll  angle 

Z  =  distance  along  trajectory 
A  and  B  =  constants  (boundary  condition) 

S  and  C  ^  =  constants  (model  construction  and  aerodynamics) 


The  values  of  these  constants  and  their  Probable  Error  as  obtained  from  the  reduction  procedure 
are  given  in  Table  EL  The  "fit"  of  the  roll  equation  (equation  1)  to  the  range  data  is  given  as  a  Probable 
Error  in  Table  IV.  Although  the  'fit"  is  very  good  since  the  estimated  error  in  angle  measurements  is 
.7°,  the  values  of  the  aerodynamic  derivatives  for  models  5  and  6,  as  obtained  from  the  constants  by 


2IC2 


(2) 

(3) 


(4) 


are  of  relatively  poor  quality  compared  to  model  7,  as  indicated  by  their  Probable  Errors  (See  Table  V). 
The  results  on  an  individual  round  depend  on  the  functioning  of  the  range.  A  small  number  and/or  a  poor 
distribution  of  working  stations  will  give  data  which  although  having  an  excellent  ’fit"  may  3deld,  as  in  this 
case,  relatively  poor  aerodynamic  derivatives. 

^BRL  Report  No.  711 
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These  new  data  are  plotted  in  figures  17,  19,  and  20  together  with  (1)  the  curve  obtained  from  linear¬ 
ized  theory  as  applying  to  this  configuration  and  (2)  the  previous  data  on  this  configuration.  In  fairing  a 
curve  through  these  data  an  attempt  was  made  to  take  into  account  the  relative  weights  as  obtained  from 
the  reciprocal  of  the  Probable  Error,  which  indicate  that  model  5  and  6  should  be  weighted  approximately 
1/4  that  of  model  7. 

The  results  show  that  the  experimental  curve  for  C  at  M  =  1.7  deviates  from  linearized  theory  by 
16%  and  that  the  deviation  increases  with  increasing  Mach  number  until  at  M  =  2.7  the  value  is  22%.  The 

experimental  curve  for  C5  at  M  =  1.7  deviates  by  7%  and  at  M  =  2.7  by  17%.  The  experimental  curve  for 

6 

6  deviates  -11%  at  M  =  1.7  approximately  -5%  at  M  =  2.7. 

One  source  for  this  deviation  maybe  the  neglection  of  the  thickness  effect  in  the  linearized  theory, 
for  if  a  comparison  is  made  to  the  lift  of  an  infinite -aspect- wing  of  such  cross  section  (16%  thickness  wing 
wedge)  computed  from  the  exact  theory  using  oblique  shock  equations,  the  result  is  an  about  15%  lower 
value  for  the  linearized  lift  at  M  =  1.7  and  about  22%  at  M  =  2.6.  Although  this  would  tend  to  explain  the 

variation  of  Cjj  it  sheds  little  light  on  which  not  only  deviates  from  theory  by  far  less  but  with  increas- 
P  ^ 

ing  Mach  number  departs  more  rapidly  from  theory  than  does  C.  .  A  brief  discussion  of  these  derivatives 

1  P 

is  given  in  the  previous  paper  but  there  is  still  no  satisfactory  explanation  of  the  variation  of  these  deriva¬ 
tives  with  Mach  number. 

TABLE  I 

Thick-Finned  Roll  Models 


Round 

Missile 

Av.  Angle  of 
Cant  in  Two 

Moment  of  Inertia 
in  Roll  „ 

Mass  of 
Missile 

No. 

No. 

Tail  Fins 

Material 

I  (slugs  -  ft  ) 

M  (slugs) 

2125 

5 

1.9333'’ 

Aluminum  fuselage 

1.266  X  10'^ 

.012256 

2324 

6 

2. 1000° 

&  bronze  fins 

11 

1.263  X  10’^ 

.012225 

2322 

7 

4.4542° 

n 

1.243  X  10“^ 

.012208 

TABLE  II 

Thick-Finned  Roll  Models 

Vel.  of  Sound 

Missile  Mach.  No.  at  Mach  No.  at  Mach  No.  at  in  Range  Air  Density  of 


No. 

Beginning  of 
Trajectory 

End  of 
Trajectory 

Center  of 
Data 

Atmosphere 

ft/sec 

Range  Atmosphere 
slugs/ft^ 

5 

2.722 

2.587 

2.682 

1125 

2.342  X  10"^ 

6 

2.301 

2.166 

2.247 

1132 

2.309  X  10'2 

7 

2.450 

2.327 

2.416 

1132 

2.311  X  lO'S 
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TABLE  in 

Thick-Firmed  Roll  Models 

Model  5  6  7 

P.E.  .7"  .3“  .6° 


TABLE  17 

Thick-Firmed  Roll  Models 
(yaw  cards) 


Model 

5 

6 

7 

A 

918.12° 

924.91 

1930.38 

B 

-919.60° 

-900.19 

-1776.07 

"^1 

.012383 

.012692 

.012786 

%  Error 

1.50 

1.98 

.36 

s 

14. 173 

15.777 

33.579 

%  Error 

.42 

.87 

.09 

c„ 

. 175504 

.200242 

.429341 

%  Error 

1.56 

2. 16 

.37 

1.836  X  10'^ 

2. 160  X  lO""^ 

1.864  X  10”'^ 

%  Error 

1.15 

.14 

.65 

^  2 

Yaw 

1.6880 

1.5894 

.0919 

14  14  19 


Number  of 
Observations 
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TABLE  V 

Thick-Finned  Roll  Models 


Model 

5 

6 

7 

Mach 

2.682 

2.247 

2.416 

(Theor.)  C 

*P 

.3100 

.3622 

.3399 

(Ex.)  Cj 

P 

.4014 

.4258 

.4242 

P.E.  % 

1.48% 

1.95% 

.35% 

Dev.  from  theory 

29.5% 

17.6% 

24.8% 

(Theor.)  C. 

.4605 

.5494 

.5107 

(Ex.)  C 

5 

.5695 

.6177 

.6192 

P.E, 

1.56% 

2.16% 

.37% 

Dev.  from  theory 

23.7% 

12.4% 

21.2% 

(Theor.)  ^ 

.02506 

.02780 

.05839 

(Ex.)  ^ 

.0236 

.0261 

.0559 

P.E. 

.43% 

.88% 

.09% 

Dev.  from  theory 

-5.83% 

-6.11% 

-4.26% 
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